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Abstract. Thermalization in highly excited quantum many-body system does not necessa- 
rily mean a complete memory loss of the way the system was formed. This effect may 
pave a way for a quantum computing, with a large number of qubits n ~ 100-1000, far 
beyond the quantum chaos border. One of the manifestations of such a thermalized non- 
equilibrated matter is revealed by a strong asymmetry around 90° cm. of evaporating proton 
yield in the Bi(7,p) photonuclear reaction. The effect is described in terms of anomalously 
slow cross symmetry phase relaxation in highly excited quantum many-body systems with 
exponentially large Hilbert space dimensions. In the above reaction this phase relaxation is 
about eight orders of magnitude slower than energy relaxation (thermalization). 
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1 Introduction 

Independent-particle modes in interacting many-body systems result from a mean-field appro- 
ximation and are at the center of many theoretical considerations. In the case of a quantum 
information device this independent-particle basis can be considered as the basis spanned by the 
individual qubits, the "computational basis". At high excitation energy, the interaction between 
the particles results in a rapid mixing of the independent particle states PJE1IH1II]- This mixing 
leads to the formation of complicated many-body configurations. Each of these individually 
ergodic (independent of the initial conditions) many-body states is characterized by sharing the 
energy between many particles of the system. The characteristic time for the formation of such 
thermalized many-body states is given by the inverse spreading width, r t h = h/r spi [HI 12- The 
quantity T spr also characterizes the width of the distribution of the expansion coefficients of the 
many-body eigenstates over a noninteracting mean-field basis p]|21IHllI]. 

There is another way to interpret spreading width if there exists a classical analog to the 
system. Then it is given by the width of the energy distribution of trajectories determined by 
the Hamiltonian with interaction, measured in terms of the independent particle Hamiltonian. 
This distribution also provides a description of semiclassical wave functions [2] around which 
the actual wave functions have Gaussian fluctuations 
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The question now arises whether phase relations and/or correlations between these indivi- 
dually ergodic, spatially extended, many-body states in the superposition may still preserve a 
memory of the way the system was excited. This question is of fundamental importance for 
the study of relaxation phenomena in nuclear, atomic, molecular and mesoscopic many-body 
systems, and for many-qubit quantum computation. In particular, we recently proposed [7j 
that, if phase relaxation is slower than energy relaxation, this can extend the time for quantum 
computing beyond the so called "quantum-chaos border" [SJ |§1 ■ 

To answer this question from first principles is presently not possible due to computational 
limitations. Indeed, in order to solve the full quantum many-body problem one would require 
a many-qubit quantum computer. Therefore, the only currently available resort to search for 
possible manifestations of long phase relaxation is the experiment, and a careful data analysis. 
Nuclear systems are an ideal testing ground to study many-body systems, since nuclear interac- 
tions are so strong that external perturbations can often be neglected. In particular, the analysis 
in Ei] of the data displaying forward peaking in angular distributions of evaporation protons 
from heavy nuclei in nucleon-induced reactions indicates that phase relaxation times can be up 
to five orders of magnitude longer than energy relaxation times. 

In this paper we analyze another reaction that indicates the formation of thermalized non- 
equilibrated matter. This new form of matter was introduced by one of us I12j. Again this 
is revealed by a strong asymmetry around 90° cm. of the evaporating proton yield, but now 
in the Bi(7,p) photonuclear reaction The experiment indicates that in this case the effect 
is even more pronounced. We shall see that here the phase relaxation is about eight orders of 
magnitude slower than thermalization, making the possibility of quantum computation beyond 
the quantum chaos border an even more attractive concept. 

The present work is a step toward a more realistic situation because the entrance channel is 
given by electro-magnetic interaction. This is more similar to loading mechanisms in proposed 
quantum computers. 

The present work does not suggest to use photonuclear reactions as a practical experimental 
setup for universal quantum computation. Indeed, one can not perform a universal set of gates 
for the considered photonuclear reaction. However the proposed analysis does demonstrate that 
phase relaxation can be much longer than thermalization which, in turn, illustrates the main 
idea of quantum computing far beyond the quantum chaos border. 

2 Experimental manifestation of thermalized non-equilibrated 
matter in Bi(7,p) evaporation process 

We analyze the proton yield of the Bi(7,p) photonuclear reaction produced by 24-MeV brems- 
strahlung. In Fig. 1 we present an angle-integrated photo-proton spectrum scaled with the 
outgoing proton energy e times the cross section cri nv (e) of the inverse process of the capture 
of the proton with energy e by the residual nucleus. This e-dependent inverse cross section is 
determined by the penetrability of the Coulomb and centrifugal barriers and was taken from |14j 
for a nuclear radius of 1.5 x lO -13 ^ 1 / 3 cm, where A is the nuclear mass number. 

One can see that, for s < 8 MeV, the scaled spectrum has an exponential shape with a slope 
of 0.55 MeV. This is characteristic for the decay of thermalized compound nucleus with a "tem- 
perature" T = 0.55 MeV of the residual nucleus. We evaluate the average excitation energy of 
the compound nucleus to be E* = 14 MeV, i.e. slightly above the center of the giant resonance 
peak at 13.5 MeV For e = 4 MeV, the average excitation energy of the residual nucleus 
is E* es = 6.3 MeV. The average number of excitons (particles plus holes) hi the residual 
nucleus is n = (2gE* cs ) 1 ^ 2 ~ 14, where g = A/ 13 MeV -1 is the level density of independent 
particle states. The standard deviation of n is about (n/2) 1 / 2 ~ 2.6. Therefore, the energy per 
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Figure 2. Experimental proton angular distri- 
butions (in arbitrary units) from the Bi(7,p) pho- 
tonuclcar reaction for e = 2-8 MeV The 
curve is a fit to the experimental data (see text). 
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Figure 1. Scaled experimental angle-integrated 
spectrum (in arbitrary units) of protons of 
the Bi(7,p) photonuclear reaction produced by 
24 MeV bremsstrahlung [Ej. The line is expo- 
nential fit of the scaled spectrum for e < 8 MeV 
with the slope (nuclear "temperature") of 0.55 
MeV (see text). 

exciton in the residual nucleus, i.e. the nuclear "temperature", is T = 0.38-0.55 MeV. This is 
consistent with the fit in Fig. 1 for e < 8 MeV with T = 0.55 MeV. Similarly we find that the 
temperature of the Bi compound nucleus with the excitation energy 14 MeV is 0.57-0.78 MeV 
which is roughly consistent with the fit in Fig. 1. This implies thermalization of both compound 
and residual nuclei |14j . 

In Fig. 2 we present experimental proton angular distributions from the Bi(7,p) photonuclear 
reaction for e = 2-8 MeV |13j . One can see that, in spite of a complete energy relaxation in the 
thermalized compound nucleus, the angular distributions are strongly asymmetric about 90°, 
i.e. memory of the direction of the incident 7-ray beam is clearly retained. 



3 Determination of the cross symmetry phase relaxation width 

In order to identify anomalously slow phase relaxation behind the forward peaking of the eva- 
porating protons we make use of a standard formula for double differential cross section in 7-ray 
induced reactions [T7] : 

a(9) = (A 2 /8vr)(2/ + l)' 1 £ S^^ErS^^E) 

(Li 1 - 1 1\L X 1 h 0)(£ 2 1 - 1 1\L 2 1 h 0)A(h,p)A(l 2 ,p) 

[(2Ji + 1)(2J 2 + l)] 1 / 2 i L2 - Ll+/l - i2 (-l) s '- / - i+Ll - i2+1 VF(JiL 1 J 2 L 2 ;/L) 

Z(l x L x l 2 L 2 - 1 L)Z{l' v hl' 2 J 2 - s' L)P L (9), (1) 

where the sum is over Li, L 2 , h, l 2 , Ji, J 2 , n±, tt 2 , l[, 1' 2 , s', p and L. Here A is the wave length of 
the electromagnetic radiation, E the total energy of the system, / the spin of the target nucleus, 
and a, a' are microstates of the reaction partners in the entrance and exit channels, accordingly. 
Ji, J 2 are the total spins of the compound nucleus, tt\, ir 2 its parities, L\, L 2 the total angular 
momenta of multipoles, l\, l 2 the orbital momenta in the entrance and l[, 1' 2 in the exit channel, 
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and s' is the channel spin in the exit channel. The symbol p is defined to have the value 
for magnetic radiation and 1 for electric radiation. For p = 0, A(/ 1 ( 2 ),p) = ^/i( 2 ),Li( 2 ) an d> f° r 
p = 1, A(/ 1 ( 2 ),p) = ^i {2 ),ii(2)±i- I n equation (-|-) denote the Clebsch-Gordan coefficients, 
W are the Racah coefficients defined in ^S], and the Z coefficients are defined in |19| . S-matrix 
elements are in a mixed representation in which the entrance channel states are in the multipole 
representation, and the exit channel states are in the channel spin representation. 

In what follows we use a symplified form of equation We neglect the proton spin in the 
exit channel so that the spin of the residual nucleus, is equal to the exit channel spin. We 
neglect the target spin in the entrance channel, 1 = 0. Therefore, J±(2) — -^1(2)- We take into 
account the proton orbital angular momenta l[ 2 < 2 only since, for ^ 2 — 3, evaporating protons 
are significantly sub-barrier due to the centrifugal and Coulomb barriers. We consider electric 
dipole (Li,2 = 1) and quadrupole (£1,2 = 2) radiation only. As a result, equation © takes the 
form 

aid) = {>?/^)Y, s X^s E y s XM E ) 

(Li 1 - 1 l|Li 1 h 0)(L 2 1 - 1 1\L 2 1 h 0)6 ll)Ll±1 6 l2iL2±1 

^-L.+h-h {-if^zihUhL^ 1 L)Z(l , 1 L 1 l , 2 L 2 ;l' L)P L (9), (2) 

where the sum is taken over L\, L 2 , 1%, Z 2 , l[, 1' 2 , I' and L. In equation ((2j), the states with 
different total spins, L\ 7^ L 2 , correspond to opposite parities. For L\ = L 2 (L\ ^ L 2 ), the sum 
in equation (j2j) includes the terms with even (odd) values of \l[ — l' 2 \ only. 

The key element in a description of the asymmetry of the evaporating protons is a correlation 
between fluctuating S'-matrix elements with different total spins L\ 7^ L 2 (L\ = 1 and L 2 = 2 
or L\ = 2 and L 2 = 1) and parities [111 I12j : 

l \ S aL 1 ; ( x'V^ E y S a L 2 ; ( x'V^ E )) ~ (l + f3/T cn ) ' ^ 

where the brackets (•) stand for ensemble averaging which is equivalent to the energy [E) ave- 
raging under the ergodicity condition. The physical meaning of the spin and parity off-diagonal 
("cross symmetry") phase relaxation width (3, introduced in |2U1 1111 IT2*] . is the characteristic 
inverse time over which the interference between the states with different total spins and parities 
does not vanish upon the energy averaging and therefore the memory of the direction of the 
initial beam is preserved [111 I12j . If this phase memory time t p ^ = hj (3 is about as long or 
longer than the average life time h/T cn of the compound nucleus then the evaporating protons 
are emitted asymmetrically about 90° cm., i.e. the memory about the direction of the initial 
beam is retained. However, if the phase memory time is much shorter than the average life time 
of the compound nucleus then the spin and parity off-diagonal correlations vanish, memory on 
the direction of the initial beam is lost leading to the conventional Bohr picture of compound 
nucleus with the symmetric about 90° cm. angular distributions. Therefore, deviation of the 
angular distributions of the evaporating particles from the symmetry around 90° cm. indicates 
that (3 is smaller or comparable to r cn . 

For the orbital momentum off-diagonal correlation between S'-matrix elements with the same 
total spins and parities we have |11| I12j 

(SL;a % m*S^ L;a% {E)) = [(IS^^I^dS^, {E)\^l\ (4) 

Both the equations © and @ reflect a strong correlation between the partial width ampli- 
tudes and, as a result, between S-matrix elements with different orbital momenta (l[ 7^ 1' 2 ) of 
evaporating protons referred to ^] as the continuum correlation. Note that such a strong 
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correlation between reduced width amplitudes corresponding to the same total spin and parity 
values but different orbital momenta was experimentally revealed for a number of compound 
nuclei in the regime of isolated resonances |21| . 
The statistical model |14j yields 



Here, = T L are the entrance channel transmission coefficients for the formation of the 
compound nucleus with the total spins L = 1 and L = 2 due to the absorption of electric dipole 
and quadrupole radiation, accordingly. The exit channel transmission coefficients are assumed 
to be independent of the compound nucleus spin L and the spin of the residual nucleus I' [Tlj . 



We use equations ©, ©, @ arid (J5J) for the analysis the experimental angular distribu- 
tions in Fig. 2. We find that, upon energy averaging and summation over microstates of the 
residual nucleus, apart from the overall normalization constant, a shape of the angular distri- 
butions depends on the four parameters: A = t l=2 /T l=1 , B = T l ' =1 /T l ' =0 , C = T l '= 2 /T l ' =0 , 
and 0/T cn . 

From the fit of the experimental angular distributions in Fig. 2 we obtain: A = 0.082, 
B = 0.47, C = 0.37 and (3/T cn = 0.11. The compound nucleus' decay width r cn for Bi with 
an excitation energy of 14 MeV can be estimated from the systematics in Fig. 7 of [22], which 
provides a good description of the experimentally determined r cn for a wide range of mass 
numbers. From this estimation we obtain r cn ~ 0.1 eV what yields in turn ~ 0.01 eV. At 
the same time, the standard nuclear physics estimate for the spreading width of Bi nucleus with 
the excitation energy 14 MeV is about 2 MeV (see Fig. 2.1 in US])- This is close to another 
estimate of T width of a dipole giant resonance 0, which is about 4.5 MeV for Bi |15j . 

It corresponds to an exponentially large effective dimension of Hilbert space given by iV e g ~ 
r sp r/£> ^ 10 16 , where D ~ 10~ 16 MeV is the average level spacing of the Bi compound nucleus 
with the excitation energy 14 MeV. Note that the total spin and parity off-diagonal S-matrix 
correlations for evaporation processes were justified in in the limit N c q — > oo. For the cross 
symmetry phase relaxation time much longer than the energy relaxation time, the formalism also 
leads to (i) quantum-classical transition |24U25| . stable coherent rotation |26U27ll25] . Schrodinger 
cat states |291 130 j in complex collisions in the regime of strongly overlapping resonances of the 
intermediate system, and (ii) spontaneous correlations, non-equilibrium phase transitions |31| 
and anomalous sensitivity in finite highly excited many-body systems |32j . 

It should be noted that while we have been able to determine an anomalously small value 
of f3 from the data analysis, its theoretical evaluation is currently an open problem. Therefore 
more theoretical insight is needed for a deeper understanding of the effect of anomalously slow 
cross symmetry phase relaxation in highly excited quantum many-body systems. 

4 An illustration of the idea of quantum computing far 
beyond the quantum chaos boarder 

The idea of quantum computing on a time scale r p h = h//3 — 6 x 10 -14 sec, which is much longer 
than the thermalization time h/T spr ~ 3 x 10" 22 sec, can be illustrated as follows [7j. Consider 
the Bi(7,p) compound nucleus reaction to be the quantum protocol. The single-particle basis is 
the quantum register. The entrance channel represents the loading process. The measurement 
of the angular distribution of the yield of evaporating protons plays the role of the readout. In 
accordance with the standard criterion jHUH], a quantum computer melts down at t > h/T spT . 
If so, then any information about specific features of the loaded state (entrance channel) is lost 
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while the protocol is executed. Yet, as was shown in the previous section, our "read out" (the 
angular distribution of evaporation protons) shows that our "protocol" yields non-trivial results 
eight orders of magnitude later than the time scale h/T spr ~ 3 x 10~ 22 sec. for the onset of 
quantum chaos in our "quantum computer". These results depend on the "loading", i.e. the 
direction of the incident beam. 

In order to demonstrate that the obtained phase memory time r p h = Kj ' j3 ~ 6 x 10~ 14 sec. 
is indeed a very long time scale of the system we note that during this time a nucleon at the 
Fermi energy crosses the Bi compound nucleus about 10 8 times. Yet, this phase memory time is 
still about eight orders of magnitude shorter than the Heisenberg time h/D ~ 10~ 5 sec, where 
D ~ 10 -16 MeV is the average level spacing of the Bi compound nucleus at an excitation energy 
of 14 MeV. 

Clearly, the spectrum of the Bi compound nucleus in the above considered reaction is not 
resolved since Y cn /D ~ 10 , i.e. resonances of the compound nucleus are strongly overlapping. 
Therefore, all exponentially large information hidden in the Hilbert space is not available. Ho- 
wever we have demonstrated that useful information, such as the phase relaxation time, still 
can be obtained even though the quantum protocol time is about nine orders of magnitude 
shorter than the Heisenberg time for our "quantum computer" . This is in accord with a similar 
observation in [§]■ 

5 Conclusions 

We have demonstrated that thermalization in highly excited quantum many-body system does 
not necessarily mean a complete loss of memory of the way the system was formed. Manifestation 
of such a thermalized non-equilibrated matter has been revealed from the analysis of a strong 
asymmetry around 90° cm. of the evaporating proton yield in the Bi(7,p) photonuclear reaction. 
We have shown that thermalized non-equilibrated matter can exist for time spans of ~ 6 x 
10~ sec, which is eight orders of magnitude longer than thermalization time in this example. 
This indicates that long lived transient states can exist in many-body systems with exponentially 
large dimensions of Hilbert space. If a quantum computer with a large number of qubits n ~ 100- 
1000 can be brought into such a state, this may provide a solution for the scaling problem which 
is one of the central challenges of quantum information |33j . 
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